Disordered graphene and boron nitride in a microwave tight-binding analog 
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Experiments on hexagonal graphene-like structures using microwave measuring techniques are 
presented. The lowest transverse-electric resonance of coupled dielectric disks sandwiched between 
two metallic plates establishes a tight-binding configuration. The nearest-neighbor coupling approx- 
imation is investigated in systems with few disks. Taking advantage of the high flexibility of the 
disks positions, consequences of the disorder introduced in the graphene lattice on the Dirac points 
are investigated. Using two different types of disks, a boron-nitride-like structure (a hexagonal 
lattice with a two-atom basis) is implemented, showing the appearance of a band gap. 



PACS numbers: 42.70.Qs, 73.22.-f, 71.20.-b, 03.65.Nk 

I. INTRODUCTION 

Graphene, due to its unusual and astonishing elec- 
tronic and physical properties, is a subject of intense 
focus, especially since it was experimentally realized in 
2005 by K. S. Novoselov, A. K. Geim and coworkers 1 , 
who were honored by the Nobel prize in physics in 2010. 
One of graphene's amazing features is the linear disper- 
sion relation close to the i^T-points and the Dirac like 
Hamiltonian whose consequences can be seen in the den- 
sity of states (for details on graphene see the Revs. H 
and0). Graphene is one of the candidates to take over 
from silicon-based electronics. Graphene has a macro- 
scopic ballistic carrier transport due to its high charge 
carrier mobility and a high magnetic susceptibility /jA. A 
problem is still that graphene in its pure form is metal- 
lic and has no band gap AE. Nevertheless, a band gap 
can be induced by using bilayer graphene^, disorder^!, 
lateral-superlattice, epitaxially grown graphene on top of 
a crystal-like boron nitride^— or by confincmcnt o 1 1 ' 12 . 
Applying stress to graphene can induce effects on the 
electron transport like a magnetic field 13 . Thus graphene 
offers a big variety of interesting possibilities from both 
applied and fundamental points of view. 

In real single layer graphene, it is quite difficult to pre- 
pare special shapes, induce precise structures of vacan- 
cies and so on. This can be easily performed in experi- 
ments with classical waves. For example in microwave 
analog experiments, the Dirac poin1jii~— , a topologi- 
cal phase transition^, the global band structures^, and 
edge state s 16 ' 18 have been observed. These states, which 
are located only at the zig-zag edges can lead to self- 
guiding unidirectional electromagnetic transport if the 
hexagonal lattice is created by ferritic scatterers^. Edge 
states have also been found in acoustics 2 ^. 

In this paper, we mimic the established tight-binding 
form of graphene's Hamiltoniart2i~— with a microwave 
setup using dielectric disks between two metallic plates. 
In a previous publication, we investigated the Dirac point 
in the density of states (DoS)^ 6 -. Here, we concentrate 
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FIG. 1. (Color online) Graphene as a superposition of two 
triangular sublattices A and B with lattice vectors a\ , d*2 and 
the vectors b%, 62, &3 connecting nearest neighbors. 

more on the experimental details of the realization of the 
tight-binding Hamiltonian (see Sec. IIIip . We show the 
flexibility of the experimental ansatz by investigating two 
different cases. The first case is to introduce positional 
disorder, where for strong disorder the Dirac point (see 
Sec. |V| vanishes. The second case is to realize a boron- 
nitride-like structure by using two different kinds of disks 
on the lattice (see Sec. I VI I) . In this case, a band gap 
opens. 



II. THEORETICAL BACKGROUND 

Graphene is a two-dimensional crystal of carbon atoms 
ordered in a regular hexagonal lattice, a so-called "hon- 
eycomb" structure (see Fig. [1]). It was first completely 
described theoretically by Wallace in 1947—. The far- 
reaching equivalence of graphene's Hamiltonian and the 
Dirac Hamiltonian for a free particle is depicted in 
Ref. [24|. Using the order and the notation of Ref. [24], 
we first give a short repetition of a standard derivation 
of graphene's Hamiltonian to underline its equivalence to 
our experiment. 

Graphene's characteristic honeycomb lattice structure 
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is caused by the sp 2 bond of the carbon atoms resulting 
in three cr-bonds in the plane including an angle of 120° 
and a free p z orbital perpendicular to the plane. The 
bonding length between the carbon atom is a& = 1.42 A. 
The basis vectors of the hexagonal lattice — illustrated 
in Fig. Q]— are ai = (VS/2, -1/2) a and a 2 = (0, 1) a, 
where a = Vc^b is the lattice constant. The vectors 
connecting the two sublattices are b\ = (1/2-^/3, 1/2) a, 
b 2 = (-1/^3,0) a, and b 3 = (1/2^3,-1/2) a. In the 
tight-binding approximation, we assume a basis of local- 
ized eigenf unction at the lattice sites imitating the p z - 
orbitals. Thus we can write down the Hamiltonian H in 
second quantization which includes the position energy 
and the hopping term. 

For the position term H pos , we use a diatomic system 
(indicated in Fig. Q] by A and B) placed on the two trian- 
gular lattices A and B, which is realized, e.g., in boron 
nitride. In the tight-binding model an atom of sort A only 
interacts with its three nearest neighbors that all belong 
to the lattice B (see Fig. [T]). Second nearest- neighbor 
hopping is here neglected. The energy difference between 
electrons localized on lattice A and electrons on lattice 
B is called j3 — (e^ — es)/2. The associated part of the 
Hamiltonian can be written as 



e A n(A) + e B n(B) 



fi^[n{A)-n{A + b 1 )] 



e~A + e B 



Ntot (1) 



where n(A) and n(A + b\) are the counting operators 
of electrons at positions A belonging to lattice A and 
at B = A + bi to lattice B, respectively. They can be 
written using the annihilation (U, V) and creation oper- 
ators (W, yt) as n(A) = W(A)U(A) and n(A + b x ) = 

V^(A + bi)V(A + bi). J2a 1S tne sum over au atoms of 
lattice A. The additional term proportional to Eq times 
the total number of electrons is only a constant and can 
be neglected in the further treatment. It corresponds to 
an energy shift to the mean energy of the two atoms (for 
graphene to the eigenenergy of one atom). 

The second part of the Hamiltonian includes the hop- 
ping of electrons from one lattice site to one of the 
neighboring sites. The hopping parameter n describes 
the probability for this process and corresponds to the 
nearest-neighbor coupling, which is accessible in the mea- 
surements (see Sec. IIIIE|) : 

H hop = K J2 [tf(A)V(A + bi) + Vi(A+bi)U(A)] (2) 

A,i 

with i = 1, 2, 3 for the three nearest neighbors in lattice B 
of an atom A. The spin of the electron is here neglected, 
as a spin is not existing in the microwave experiment. 
The entire Hamiltonian equals the sum of Eqs. ([T]) and 
@ and has the eigenenergies: 



E(k) = ±\/3 2 + k 



oik-b-2 



For the case of graphene (f3 = 0), the two bands with plus 
and minus signs are touching each other at the roots of 
the sum of the three exponential functions. These are 
given at K — J? - (1, a/3) and K' = —K, the corners 
of the Brillouin zone. The dispersion relation is approx- 
imately linear close to these K points with the so-called 
Fermi velocity Vf as the proportionality constant: 



E(K + q) 



for \q\ < \K\ (4) 



This equals the relativistic dispersion relation E(p) = 
±c|p| for a massless particle, when Vf plays the role of 
the speed of light, and leads to the common name "Dirac 
points" for K and K' . 

The strong connection to the Dirac theory is, moreover, 
given by the equivalence of the Dirac Hamiltonian and 
a certain form of graphene's Hamiltonian, which can be 
derived via a basis transformation described in Ref. [24]. 

From Eq. the density of states (DoS) in the vicinity 
of the Dirac points can be calculated: 



P(E) 



V3n k 2 ' 



(5) 



which is again linear. 

In our experiments, the DoS (strictly speaking the local 
density of states, LDoS) is an accessible quantity. In our 
previous work, we have observed the linear dependence of 
the DoS close to the Dirac points. In this paper we want 
to investigate the stability of the Dirac points against 
disorder (see Sec. Iv|). 

If the on-site energies arc different, i.e., (3^0 like in 
boron nitride, then a gap is generated at the Dirac point, 
where the gap width is given by 2/3 = — e B ■ An analog 
to boron nitride is presented in Sec. IVI1 



(3) 



III. EXPERIMENTAL REALISATION OF A 
TIGHT-BINDING SYSTEM 

A. Introducing microwave experiments 

Microwave experiments have been proven to be an 
attractive implementation to perform quantum analog 
measurements in the realm of "quantum chaos" for 
20 year o 25 i 26 . Exploiting the formal analogy of the 
Helmholtz and the Schrodinger equations, they have been 
used to investigate many different effects like eigenvalue 
distributions, nodal domain s 27 ' 28 , affirmation of peri- 
odic orbit theor y 29 ' 30 , random matrix theor3 r 31 i 32 , and 
the random plane wave model 3 ^. Quantum-mechanical 
currents and the visualization of chaotic wave func- 
tions in billiard system s 34 ' 35 have already been in the 
center of attention. Nowadays, the area of applica- 
tion became much wider and reaches beyond the first 
investigations on chaotic billiards. For instance, one- 
dimensional wave guides with introduced disorder have 
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been suitable to check scattering theory phenomena. Ef- 
fects of uncorrected and correlated disorder have been 
observed^— and localization phenomena have been in- 
vestigated in two-dimensional systems using dielectric 
disks^. Photonic band gap materials (PBGM) (for 
a review, see, e.g., Ref. 1401 ) assuming infinitely long 
cylinders have been suggested for two-dimensional tight- 
binding realizations 41 . Effects of disorder have been in- 
vestigated theoretically 4 ^ and implemented successfully 
experimentally 4 ^. 

Transport of electrons in lattices is related to wave 
phenomena. As large interest in the special properties 
of the band structure of graphene occurred, microwave 
analog experiments have been performed on graphene- 
like structures using different kinds of realizations 1 ^ -16 ' 18 . 
Here, we want to emphasize the possibilities of using di- 
electric disks with a high index of refraction^ 6 ". This new 
field offers a high flexibility to realize random or symmet- 
ric lattice structures and investigate its characteristics by 
observing the local density of states at different positions. 




FIG. 2. (Color online) Spectrum of a single disk: black line 
corresponds to a measurement with the TE antenna shown in 
the inset; orange line shows the same measurement performed 
with a TM antenna. The arrow marks the first TE resonance 
that will be used for further experiments (from Ref. [l6h . 



B. Setup 

For the investigation of tight-binding systems we need 

wave functions that exists in a finite volume and decrease 

faster than t4o outside. This scenario can be consid- 
er _ 

ered as a potential well configuration. For the realization 
we chose dielectric disks (Temex-ceramics, E2000 series) 
with a dielectric permittivity e r 36 and a relative per- 
meability /i,. = 1 . They have a diameter of 8 mm and a 
height of 5 mm (see inset of Fig. [2]). 

These disks were previously used to realize disordered 
systems^. They are placed on a flat copper plate. Atop 
another metallic plate is positioned parallel to the bot- 
tom plate. The two metallic plates are separated by a 
distance h. Thus the cutoff frequency for the first TE 
mode in air due to the plates is given by v cu t = co/(2/i), 
where cq is the speed of light in air—. The investigated 
resonance frequency of the disks is smaller than i/ cut guar- 
antcing that the TE waves are evanescent outside the 
disk. This assures additionally the two-dimensionality of 
the system 26 . Further descriptions about the disk eigen- 
modes are given in Sec. IIII CI 

In Fig. [21 the reflection spectra are shown using two 
different kinds of antennas. The first is a pure dipole 
antenna, which consists of a vertical part only, exciting 
only TM modes (shown in orange). The second has an 
additional horizontal part that was soldered to the end of 
a standard dipole antenna (see inset) and can also excite 
the TE modes (spectrum shown in black). The sharp 
resonance marked by the arrow is the TEi resonance of 
the disk as it is not observable with a dipole antenna. 
Besides, it is nicely separated from the neighboring res- 
onances and has a quite flat background. Thus it is a 
suitable resonance for further investigations. We always 
make sure that the other resonances are still separated 
from the frequency range of this TE resonance so that 



mixing or overlapping of TE and TM modes are avoided. 
Consequently, we work in a one-level regime, where each 
disk brings in only a single resonance. 

To mimic molecules or solids several disks have to be 
placed. As we want to mimic molecules consisting of one 
type of atoms, we have to ensure that each disk has the 
same eigenfrequency. One problem for this realization 
is the sensitivity of the resonance frequency once the top 
plate is very close to the top of the disks. This can be seen 
in Fig. [3j where the resonance frequency and width are 
shown in dependence of the height h t . ht is the distance 
between the top of the disk to the metallic top plate. At 
ht=0 mm, the bottom plate is pressed against the upper 
plate using a step motor. Thereafter it is moved in steps 
of 0.1mm. At the beginning, this just corresponds to a 
release of the pressure and thus the resonance frequency 
does not change. Otherwise, for small heights the change 
of the resonance frequency is large for even tiny variations 
of the height and than saturates for larger heights. As- 
suming height variation along the setup of a few tenth of 
a millimeter and a typical resonance width of r= 2 MHz, 
a height h t above 5 mm ensures that the variation of res- 
onance frequencies due to height variations stays within 
the resonance widths. Note that in Fig. [3] the resonance 
frequency is given in GHz whereas the width is given in 
MHz. The width decays for small distances and stabilises 
also but increases strongly once the eigenfrequency gets 
closer to the cutoff frequency. For all further measure- 
ments, a height of /i=16mm (/i t =llmm) is chosen to 
guarantee sharp and stable resonances. 

The disks are positioned by a motorized xy-stage that 
is controlled by a PC programme. Hence it is possible 
to create structures of a few hundreds of disks with an 
accuracy of the order of 0.1mm. To ensure good and 
equivalent contact of the disks with the bottom plate the 
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FIG. 3. (Color online) Resonance frequency v (upper) and 
width F (lower) of a single disk are shown in dependence of 
the distance h t between the metallic top plates and the top 
of the disk. The dashed blue line corresponds to the cutoff 
frequency at the corresponding height and the vertical dotted 
red line indicates the height at which the resonance frequency 
is above the cutoff frequency. The arrow marks the chosen 
distance ht = 11mm between the top of the disk to the top 
plate for the experiments. 



disks are pushed down while being placed. 



C. Disk eigenmodes 

The aim of our experiment is to realize a tight-binding 
system of localized functions with weak coupling to the 
outside. For the theoretical description, let us assume 
that h t =0 mm, i. e. the top plate touches the disk, the re- 
fractive index is n = 6, and the disk height is h — 5 mm. 
The TEi mode consists only of the z component of 
the magnetic held B = (0,0, B z ) and the perpendic- 
ular components of the electric field. For the given 
distance between top and bottom plates the cutoff fre- 
quency ^ cu t = 9.37 GHz for this first TE-mode. The 
z-dependence can be separated and if the resonance fre- 
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FIG. 4. (Color online) The spectra of a single disk are shown 
for different perturber positions. All spectra are obtained 
without top plate. Lowest (black) spectrum: no perturber, 
next (blue) spectrum: perturber close to the border, the fur- 
ther up the spectrum is shifted, the closer is the perturber 
placed to the disk's center. The spectra are shifted in y in 
steps of 0.05 for better visualization. The inset on the upper 
right shows a photograph of the disk including the metallic 
perturber and the measuring antenna. 



quency is below the cutoff frequency the lowest TE reso- 
nance then looks like 



B z (x,y,z,k) = B -, 



a sin (^z) K (^r); r > r D 



(6) 



where k±_ = \j k 2 - (^) , 7 = \J (f- ) - k 2 , and r = 

\/ x 2 + y 2 is the distance from the center of the disk. Jo 
and Kq are Bessel functions and a is a constant to be 
determined from the continuity equations at the surface. 
At the resonance frequency of the disk, the wave num- 
ber for r > rr> is thus purely imaginary leading to an 
evanescent wave in the area outside the disk. Note that 
in the experiment there is a gap between the top plate 
and the disks, and this description is only valid approxi- 
mately. We still use Eq. © but with an effective 7 to be 
determined from the experiment. 

To experimentally test the wave function inside the 
disk, we measured the reflection for a single disk as a 
function of the position of a perturber, i.e. a small metal- 
lic rod placed above the disk (see inset of Fig. 0} . For this 
measurement the top plate is removed. The perturber 
induces a frequency shift, which is proportional to the 
intensity of the wave at the position of the perturber. 
In Fig. 21 a monotonous resonance shift is observed con- 
firming that the wave function corresponds to the lowest 
Bessel function Jo. A direct coupling of the TEi mode 
to the TM and TMi modes is strongly suppressed by 
the continuity conditions of the fields at the disk bound- 
aries. Consequently, the wave function in the disks can 
be identified with an electron wave localized at a lattice 
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FIG. 5. (Color online) Resonance positions z/ cxp versus dis- 
tance d a between the disk and antenna (black crosses) with fit 
of Eq. (orange line) with K a =-13.89 GHz, 7=0.52 mm" 1 , 
and i^o=6.67 GHz. 



FIG. 6. (Color online) Two-disk spectra for increasing dis- 
tance. They are up-shifted with increasing distance. The 
lowest spectrum corresponds to a distance of 0.6 mm, the top- 
most to 7.6 mm. 



site. 



E. Extracting the coupling parameter in two-disk 
measurements 



D. Impact of the measuring antenna 

For all the measurements we need to know the im- 
pact of the antenna on the disk eigenmode. Hence the 
coupling between two different objects, i. e. disk and an- 
tenna, is investigated. The pure presence of the antenna 
disturbs the symmetry of the system, but the coupling 
behavior can be quantified in a simple measurement se- 
ries: we use a single disk with radius m and increase its 
distance d a to the antenna. d a is the distance to the bor- 
der of the disk and not to the center. Let us focus now on 
the shift of the resonance. The corresponding resonance 
frequencies are shown in Fig. [5] For smaller couplings, 
i. e. for larger distances, the resonance frequency becomes 
higher and its amplitude decreases. The resonance fre- 
quency approaches the pure eigenfrequency for vanishing 
coupling, which thus equals 6.67 GHz. The shift of the 
resonance due to the antenna is proportional to the wave 
function's intensity at the position of the antenna, i. e. at 
distance td + d a . Thus we fit the function 



v a (d a ) = n a \K [7 (r D + d Q )] | 2 + vq 



(7) 



to the experimentally measured resonance frequency v a . 

In Fig. [5j the experimental values are shown. A good 
agreement between fit (see orange line) and data points 
(crosses) is found. We obtained the fit parameters as 
k q =-13.89GHz, 7=0.52 mm" 1 , and ^ =6.67 GHz. It's 
worth noting that K a also includes the global coupling 
characteristics of the antenna. For an infinite distance 
the pure eigenfrequency v n = 6.67 GHz remains and de- 
fines the eigenenergy of an isolated disk. Let us now 
examine the disk-disk coupling. 



A system consisting of two disks with slightly different 
eigenfrequencies vq and v§ + A^o and a distance depend- 
ing coupling Ki(d) can be described by the Hamiltonian 



H(d) 



Kx(d) 
ki (d) p + Ai/q 



(8) 



where the chosen basis (1,0) and (0, 1) means "electron" 
in eigenstate of disk 1 with eigenfrequency isq and eigen- 
state of disk 2 with eigenfrequency isn + Avn, respectively. 
The difference Ai>q of the eigenfrequencies includes the 
influence of the antenna, quantified in Sec. lIIIDl and the 
discrepancy between the eigenfrequencies of the disks. 
The Hamiltonian has the eigenvalues 



m,i{d) = v (d) + ^-Av ± ^yji K X {d) 2 + Ai/, 



2 2 
whose difference is given by 



Av{d) = ^4 Kl (d) 2 



Ai 



The coupling parameter k,\ can be estimated from the 
overlap of the wave function of the single disks. Using 
the evanescent waves [Eq. ([6]) for r > r£>] outside the 
disks the coupling is approximately given by 



Ki(d) = k 



Kn 



70d + 2 d ) 



For the shift of the resonance we thus obtain 



Au(d) = 4/4 Kg 



Kn 



Avl 



(9) 
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FIG. 7. (Color online) Resonance splitting of the two-disk 
system depending on the disk distance (black crosses). The 
experimental data agrees well with the fit function (|10[) for 
values K = 1.35 GHz, 7 = 0.29 mm -1 , and Av = 0.021 GHz. 
The inset shows exemplary the spectrum for d — 3.3 mm 
(from Ref.Ql). 

To check it experimentally we varied the coupling be- 
tween two disks by increasing their distance in steps of 
0.2 mm and measured the splitting of the two resonances. 
Figure O shows the corresponding spectra for distances 
of 0.6 to 7.6 mm between the borders of the disks. With 
increasing distance, the decrease of the resonance split- 
ting due to weaker coupling is clearly observable. Af- 
ter fitting the resonances with Lorentzians, we can ex- 
tract the resonance splitting depending on the distances 
of the disks (see Fig. [7] crosses). The experimental find- 
ings are described by Eq. (fTTJj) and by fitting we obtained 
k =1.35GHz, 7=0.29mnr 1 , and A^ =0.021 GHz. 

In the next section we quantify the second nearest- 
neighbor coupling in a three-disk measurement. 
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FIG. 8. Coupling parameters for the three-disk measurement. 
Upper figure shows the nearest-neighbor coupling m with a 
fit of Eq. (|9]). The dashed line corresponds to the coupling pa- 
rameter ki from the two-disk measurement. The lower figure 
presents the extracted values for the next nearest-neighbor 
coupling K2- 



F. Extracting the next nearest-neighbor coupling 
parameter 

We performed an analogous measurement with three 
disks in a row whose distances will be equally increased 
in 0.2 mm steps. The antenna remains close to the cen- 
tral disk. For simplicity, we neglect in the modeling of 
such a situation the influence of the antenna and assume 
perfectly equal disks. This leads to a Hamiltonian like 

(v Kl(d) K 2 (d) 
K 2 (d) Kx(d) v Q 

with K\ the nearest-neighbor coupling and k 2 the second 
nearest-neighbor coupling. From its three eigenvalues 
"1,2,3 the expressions for the coupling can be calculated: 



Ki(d) 



V\V 2 



K2{d) = - (fi - 2 v 2 + v 3 ) 

We measured the spectra of the equidistant three-disk 
system as a function of the border to border distance d 
and extracted the corresponding eigenfrequencies. From 
those eigenfrequencies, we calculated the coupling pa- 
rameters, which are presented in Fig. [5J Additionally, 
we plot the extracted coupling coefficient k± from the 
two-disk measurement as dashed line. Overall the de- 
viations are small. Note that we have neglected in the 
extraction of «i for the three-disk system any deviations 
of the resonance frequencies Avo of the disks. Fitting the 
next nearest-neighbor coupling with the same function as 
in Sec. IIII E[ we got 



Ki{d) 
GHz 



2.93 



().:',(; j 4+i — 

2 mm 



0.009, 



3 V2 



which are similar to the fit values obtained for Ki in the 
two-disk experiment. The deviation of the fit values are 
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FIG. 9. (Color online) Ratio of coupling parameters for in- 
creasing distance of the disk. The second-nearest-neighbor 
coupling reaches only less than 15% of the nearest-neighbor 
coupling (horizontal line). 



FIG. 10. (Color online) Spectra of hexagonal cells (benzene) 
for different distances. The spectra are up-shifted with in- 
creasing distance. The lowest spectrum corresponds to a dis- 
tance of 0.5 mm, increasing in steps of 0.2 mm. 



coming from the different treatment of the differences of 
the eigenfrequencies of the disks A^. In the two-disk 
problem, they are taken into account properly, whereas 
for the three-disk system, they are incorporated in the 
three fit constants. Note that in the limit for infinite 
distances, the coupling does not go to zero but to a con- 
stant value. The agreement of the coupling constants in 
case of the two- and three-disk systems underlines the 
applicability of the tight-binding approach. To justify 
the nearest-neighbor coupling approximation, we inves- 
tigate the ratio \k2/ki\, which is presented in Fig. [9] 
The ratio is always smaller than 15% thus the effects of 
next nearest-neighbor couplings are small but not negli- 
gible. While generally ignored in graphene, next-nearest- 
neighbor and even next-next-nearest-neighbor couplings 
are of the order of 5% and can be even larger in bilayer or 
doped graphen o 22 ' 45 ^ 46 . Including higher-neighbor cou- 
plings can, for example shift the Dirac points or generate 
asymmetric band structure a 22 i 45 ' 46 . 

In contrast to the PBGM, which are typically real- 
ized by infinitely long dielectric cylinders, the trans- 
port between the disks is evanescent, whereas in PBGM 
the waves are freely propagating between the dielectric 
cylinders 41 - - — . Thus higher-order couplings are larger in 
the PBGM as the coupling is typically decreasing with 
d~ 2 , whereas in our case, it decreases with the modified 
Bessel function \K (d)\ 2 . 



G. Benzene — the main component of graphene 

The next measurements were performed on a hexag- 
onal cell. In the following, it will be called "benzene" 
though only the 6 p z orbitals of the carbon atoms are 
imitated and their hydrogen atoms are ignored. For sym- 
metry considerations, it does not make any difference. 
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FIG. 11. (Color online) Ratios K3/K1 (crosses) and /12/fti 
(squares) vs the distance. From a distance of 2 mm, the ratios 
are smaller than 0.15 (horizontal line). 



First, we check the predictions of group theory: as ben- 
zene has the dihedral symmetry Dq^ , we expect four reso- 
nances, consisting of two singlets and two doublets. The 
highest and lowest resonances are not degenerate. We 
measured spectra of benzene for different border to bor- 
der distance d of the neighboring disks increasing in steps 
of 0.2 mm. The different spectra are shown in Fig. [TU] 
and for each we find the expected four resonances. The 
degeneracy of the two central resonances was shown in 
Ref. [TU. In case of benzene, each disk has two nearest 
neighbors two next-nearest neighbors (^2) and one 
next-next-nearest neighbors (K3). From the spectra, we 
extract the four resonances and calculated the coupling 
constants from the obtained resonances by 

1 / 

«i = g {-vo -V1+V2 + VZ), 
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K 2 = ~ ( + VQ -VI-V2+ ^3) , 

«3 = g + - 2z/ 2 + f 3 ) 



We obtain a fit function of the following form for n\: 

2 



system: 



GHz 



2.28 



0.33 4 



1 d \ 

2 mm / 



0.0068. 



The K\ for the two-, three-disk and the benzene mea- 
surements shows only small deviations (typically below 
±0.01 GHz). Again, we have a look at the ratios of 
the three coupling parameters, which are illustrated in 
Fig. [TT] Note that the ratios K2/K1 cannot be directly 
compared to the one obtained for the three-disk mea- 
surement as the next-nearest-neighbor distances are not 
the same. We observe that for distances larger than 2 mm 
both ratios are smaller than 15% (horizontal dashed line). 
As we want the influence of the higher-order neighbor 
couplings to be small, we consequently work with disk 
distances above 2 mm and for the following measurements 
on graphene flakes, we use d=4mm. 

Summarizing the basic measurements on few disks we 
conclude the following. The group theoretical predic- 
tions are fulfilled verifying the symmetry considerations. 
A consistent coupling behavior is found with small sec- 
ond and third nearest neighbor coupling. Thus we have 
realized a possibility to investigate tight-binding systems 
with a table top experiment, where several parameters 
can be varied easily. Now we proceed with hexagonal 
lattices. 

Experimental results on graphene using the same setup 
have already been published in Ref. [l|| where the Dirac 
point including the linear dispersion relation close to it 
have been observed. Additionally, edge states at zigzag 
edges and the lack of edge states at armchair edges have 
been shown. Here, we will proceed by showing the flexi- 
bility of the experiment and investigate the effect of posi- 
tional (coupling) disorder on the Dirac point. Addition- 
ally, it is possible to use different kinds of disks, where 
we will realize a correspondence to boron nitride by using 
two types of disks on the two triangular sublattices. 



IV. LOCAL DENSITY OF STATES 

As we want to relate our findings for larger systems to 
the density of states, wc derive here the relation between 
the measured reflection signal Su at an antenna to the 
LDoS. 

In quantum mechanics, the LDoS C(r, E) is defined by 



£(r,£) = ]T|0„(r)| 2 8(E - E n ) 



(11) 



where <p n (r) is the n th eigenfunction at position r and 
can be expressed via the Green function G(f, r, E) of the 



C(r,E) = -i ImflimX: 

V n 



E — E n 



lm[G(r,f,E)} 



(12) 



To relate the Green function to the reflection R 2 , let us 
start with the probability T 2 for the incident wave to be 
transmitted via the antenna into the system and escape 
the system by other means, i.e. not via this antenna. 
In our case this corresponds to the absorption mainly in 
the top and bottom plate but also within the disks. By 
energy conservation, it is connected to the reflection R 2 
and given by 



1 = R 2 



T 2 = \S 11 \ 2 



I" = l-|<Siil = <P3dS 



(13) 



where ds describes a contour around the antenna (assum- 
ing a two-dimensional setup), j is the two-dimensional 
quantum-mechanical probability current, 

j = — Im (i/)*V xy ip) 
m 

which in electrodynamics corresponds to the two- 
dimensional Poynting vector. For TE modes it is given 
by^ 



S X y 



Im (B* z V xy B z 



(14) 



We neglect the prefactors because only the proportional- 
ity of the current to lm("0* \7 xy ip) is important. Assum- 
ing a point-like source at R, the wave function if) = B z 
is proportional to the Green function G(r, R, E), as both 
are the solution of the Helmholtz differential equation 
with delta inhomogeneity. The proportionality coefficient 
A(E, R) depends on the coupling of the antenna to the 
system. Combining Eqs. (fT2"l) -Eq. (fFfj). we obtain: 



1 - \S n 



-\A(E,R)f 



xlm 



G*(r,R,E)V xy G(r,R,E) 

cx \A(E,R)\ 2 Im &(R,E) 
oc \A(E,R)\ 2 C(R,E), 

where the renormalized Green function 



&(R,E) - lim 



(15) 



removes the singularities at the antenna positions R. f3 is 
the scattering length of the antenna. A detailed deriva- 
tion of the renormalized Green function can be found in 
Ref. 113. 
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FIG. 12. (Color online) The inset on the upper right shows 
the graphene flake of hexagonal shape with maximal depo- 
sition of a— 4mm. The antenna position is marked by the 
red cross. The spectra of the hexagonal graphene flakes with 
different maximal random depositioning a are shown. The 
lowest (black) spectrum corresponds to the unperturbed one 
(cr=0mm). The spectra with cr=l, 2, 3 and 4mm are shifted 
upwards correspondingly. The bars indicate the region of av- 
eraging used in Fig. 1131 

It is not clear how the coupling parameter A(E, R) 
depends on E or even on the LDoS. Therefore it is neces- 
sary to find a proper coupling theory. With the methods 
developed in Ref. H3, it is possible to derive an approx- 
imative description of the coupling for small cylindrical 
antennas showing an almost frequency independent pref- 
actor A(E,R). Hence we finally obtain the relation be- 
tween the LDoS and the reflection that is given by 

C(R,E) oc 1 - \Sn(R,E)\ 2 (16) 

Based on this knowledge we will plot in the following 
figures 1 — |Sii| 2 on the y-axis. 

V. DISORDERED GRAPHENE 

Graphene's astonishing DoS is a consequence of its 
highly symmetric hexagonal lattice structure. But in the 
real world the perfect hexagonal symmetry of graphene 
will be distorted by different means. For example, exter- 
nal strain^, vacancies^, matching to substrates^, edges, 
and corners will destroy the symmetry. Additionally, 
they can introduce disorder, e.g., via the position of va- 
cancies in the bulk& or via the edges in finite systems^. 
Effects of Anderson localization have been numerically 
observed in transport through graphene nanoribbons 
with rough surfaces^. In a work on PBGM, it was shown 
that a band gap can be quite stable against positional dis- 
order if the original band structure comes from the Mie 
resonances^,. This corresponds to our experimental re- 
alization apart from the fact that we have two touching 
bands at the Dirac point and not a proper gap. 



One kind of perturbation that can be easily mod- 
eled by our setup is bulk disorder within the disk po- 
sition. Starting from the perfect graphene lattice posi- 
tion of the disks r^o, we create a new position fi t „ = 
ri.o + Ar"i(cos ^e x + sin (9^), where Oi is randomly cho- 
sen from {0,27r} giving the direction of the shift and Ar^ 
is randomly chosen from {0,<r} giving the length of the 
shift. e x and e y are the unit vectors in x and y directions. 
The positional disorder corresponds in the tight-binding 
description to a disorder in the coupling. 

We want to test how stable the Dirac point is against 
an introduced disorder. To reduce effects of edge and cor- 
ner states we chose a hexagonal graphene flake as shown 
in the inset of Fig. [T^l The antenna position where the 
reflection Su is measured is marked by a red cross. Start- 
ing with a regular lattice with a disk to disk distance of 
d=A mm, we add noise to every disk position as described 
above. We realized four disordered samples with a max- 
imal distance to the original position of a = 1, 2, 3, and 
4 mm. The inset of Fig. [T^J corresponds to the maximal 
value of ct=4 mm. The behavior of the LDoS is shown in 
Fig. 1121 For a = 1 mm we find a quite similar spectrum 
and can still see a strong reminiscence of the Dirac point 
in the LDoS (red arrow). For bigger shifts the structure 
disappears and the Dirac point becomes unobservable. 

To investigate this in a more quantitative manner, we 
perform averages over small frequency windows. The av- 
eraging windows are indicated by bars in Fig. [12] and are 
13 MHz at the Dirac point and 26 MHz inside the bands. 
To observe effects of the Dirac points in the averaged 
quantity, the averaging range around the Dirac point is 
chosen so small that a reduction of the DoS is still clearly 
observed. It is necessary to introduce the average as we 
are dealing with finite systems. Thus the experimental 
value at the Dirac point is dominated by the frequency 
distance of the closest resonance and its width. If small 
disorder is introduced, we are not testing the stability of 
the global structure of the Dirac point but only the sta- 
bility of the closest eigenvalue. By averaging locally, we 
wash out this effect but still keeping the information of 
the diminution of the DoS close to the Dirac point. 

In Fig. [T31 on the upper part the averaged LDoS ((1 — 
|<Sii| 2 )) is presented around the Dirac point (squares) 
and within the two bands (triangles and crosses). One 
observes with increasing disorder an increase at the Dirac 
point, whereas the averaged LDoS within the bands os- 
cillates. At a— 2 mm, the averaged values are of the same 
order, thus indicating the final vanishing of the Dirac 
point. That the Dirac point is quite stable against dis- 
order perturbations is in agreement with the findings for 
PBGMi^. 

Additionally, we performed numerical tight-binding 
simulations on the same flake as in the experiment in- 
cluding more than only nearest-neighbor couplings. We 
diagonalized the Hamiltonian 

H th = i/q 1 + K (17) 

where v§ — 6.65 GHz is the eigenfrequency of a single 
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FIG. 14. (Color online) Two spectra of the same lattice con- 
figuration, but consisting of one type of disks (graphene, or- 
ange), respectively, two types of disks (boron-nitride, black). 
The disk close to the measuring antenna was not moved, to 
guarantee a consistent coupling. The dotted lines correspond 
to the eigenfrequencies of the two types of disks and indi- 
cate the band gap size. The spectra were shifted by its corre- 
sponding Dirac frequency vo such that the center of the bands 
are the same. Additionally, the red dashed lines indicate the 
linear behavior of the DoS (see also Ref. [IrJl . The dashed 
horizontal line corresponds to the background signal of the 
experiment. 



VI. IMITATING BORON NITRIDE 



FIG. 13. (Color online) On the upper figure the frequency 
averaged local density of states (1 — 1 5"n | 2 ) around the Dirac 
point (squares) and within the bands at lower (crosses) and 
higher (triangles) frequencies than the Dirac point are shown 
for different maximal disorder a. The windows for averaging 
are indicated in Fig. [T2j On the lower part the numerically 
obtained averaged DoS is presented. Symbols are used corre- 
spondingly to the upper figure. 



disk. The coupling matrix IK contains the coupling ele- 
ments k. A disk is coupled to all neighbors that have a 
border to border distance d less than 18 mm. The cou- 
pling parameter k is calculated by Eq. (|9]) using the spe- 
cific disk distances d. The eigenfrequencies are obtained 
by diagonalizing the Hamiltonian. Then the DoS is cal- 
culated. On the lower part in Fig. [13] the averaged DoS 
for different frequency ranges are shown. The average is 
performed over the same frequency windows as in the ex- 
perimental case. In contrast to the experiment, the DoS 
is shown and not the LDoS. One observes also here that 
the density close to the Dirac point is increasing. Each 
step actually corresponds to one additional state inside 
the averaging region. Within the bands, a decrease of the 
averaged density of states is found. Overall, the Dirac 
point seems to be quite stable against the introduction 
of small white noise disorder. 



Boron nitride has a hexagonal lattice structure, but 
with two different sorts of atoms placed on the two trian- 
gular sublattices. Double layer of boron- nitride graphene 
lattices are a candidate to have on the one hand a band 
gap of similar energy than silicon and on the other hand 
keep the high charge carrier mobility of graphene^. Thus 
it is one candidate to substitute silicon and overcome es- 
pecially problems creating high frequency transistors. 

The energy difference (3 of the two types of atoms 
causes a widening of the band gap between the 
conducting and the valence band, corresponding to 
Eq. ([3j). We imitate this situation using two charges of 
disks with mean eigenfrequency of 1^=6.637 GHz and 
j/b=6.672GHz and create the same lattice once with 
equal disks and once with the different disks. A com- 
parison of the resulting LDoS is presented in Fig. [T4j 
The orange curve, already presented in Ref. [TH [see Fig 
4(a)], corresponds to a graphene lattice with the Dirac 
point at 1^=6.637 GHz. The curves were shifted by their 
corresponding Dirac frequencies (vq) so that the central 
frequency for graphene and boron nitride coincides. The 
dashed horizontal line corresponds to the experimental 
background signal, mainly coming from the coupling of 
the antenna to the TMo mode and thus corresponds to 
the for the DoS. The red dashed lines indicate the lin- 
ear behavior of the DoS close to the Dirac point. The 
black curve shows the LDoS in case of boron-nitride. A 



11 



gap and a corresponding shift of the gap are observed. 
The shift is due to the fact that we use here two differ- 
ent disks. Thus the center of the band is expected to be 
at \ (6.637 + 6.672) GHz = 6.655 GHz, which is in good 
agreement with the experiment. The gap width is ex- 
pected to be Av g = |(6.637 - 6.672) | GHz = 0.035 GHz, 
which corresponds to the difference between the dotted 
lines. It is in agreement with the experimentally observed 
gap. 

VII. CONCLUSIONS 

In this paper, we present experimental results of hexag- 
onal tight-binding configurations in a microwave setup. 
First, we show in detail how a tight-binding setup is real- 
ized. For setups with few disks, we compare next-nearest- 
neighbor coupling to the nearest-neighbor coupling. The 
flexibility of the experiment allows to vary easily disks 
positions, i.e. the coupling, but also the "on-site energies, 
i.e. eigenfrequency of the disks. We choose exemplarily to 
investigate disordered graphene. The disorder measure- 
ments show that the Dirac points seems to be quite stable 



against white noise disorder. As a second example show- 
ing the variation of the on-site energy, we observe a band 
gap by using two different type of disks corresponding to 
a boron nitride sheet. 

Exploiting the flexibility of the setup the realization of 
many different tight-binding structures are imaginable. 
Realizing different kinds of vacancies including disorder, 
transport properties along lattice axes. By lattice defor- 
mation, one can either play with K and K' point o 51 ' 52 , 
which have been experimentally realized using a similar 
setupii, or simulate pseudo-magnetic fields to generate 
associated Landau levels^. 
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